Onset of Localization in Heterogeneous Interfacial Failure 
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We study numerically the failure of an interface joining two elastic materials under load using a 
fiber bundle model connected to an elastic half space. We find that the breakdown process follows 
the equal load sharing fiber bundle model without any detectable spatial correlations between the 
positions of the failing fibers until localization sets in. The onset of localization is an instability, not 
a phase transition. Depending on the elastic constant describing the elastic half space, localization 
sets in before or after the critical load causing the interface to fail completely, is reached. There is 
a crossover between failure due to localization or failure without spatial correlations when tuning 
the elastic constant, not a phase transition. Contrary to earlier claims based on models different 
from ours, we find that a finite fraction of fibers must fail before the critical load is attained, even 
in the extreme localization regime, i.e. for very small elastic constant. We furthermore find that 
the critical load remains finite for all values of the elastic constant in the limit of an infinitely large 
system. 

PACS numbers: 81.40.Np,81.40.Pq,62.20.mm,83.80.Ab 



The joining of interfaces, e.g. by welding or gluing, is 
an important part of everyday technology, a technology 
that has been refined through the centuries. When joined 
interfaces are subject to excessive loads, failure occurs. 
Often it is not the joints themselves that fail, but the 
material that surrounds them, as the joints themselves 
are the stronger. 

The aim of this work is, however, not to study failure 
of such joints with improvement of technology in mind. 
Rather, we take the point of view that failure of a het- 
erogeneous joined interface provides a simplified model 
for fracture in bulk materials. Such an idea is not new. 
Schmittbuhl et al. Q and Schmittbuhl and Mal0y Q 
studied, first computationally, then experimentally, the 
roughness of the fracture front moving through a sin- 
tered interface between two Plexiglas plates that are be- 
ing plied apart in a mode-I fashion. The study of the fluc- 
tuations of this fracture front provides much insight into 
the much more complex morphology of three-dimensional 
fracture surfaces [3[. 

We focus on the phenomenon of localization in this 
work. Local failure occurs either because the material is 
weaker at that spot or because it is more loaded there 
than elsewhere. Differences in local strength is due to 
heterogeneities. Differences in loading is due to struc- 
ture in the stress field. If we assume that the interface is 
loaded uniformly — the local stress field will be quite uni- 
form. Local failure will occur because of material weak- 
ness. For the simple reason that the further away we 
search from a point where a failure has occurred, the 
weaker the weakest spot we have found so far will be, 
localization is disfavored. Heterogeneity in strength in- 
duces a "repulsion" between the local failures. However, 



when failed areas build up, local stress is concentrated 
at the rim of the failed areas making these regions li- 
able to fail. Heterogeneity in the stress field induces an 
"attraction" between local failures. 

Localization occurs when attraction wins over repul- 
sion. A transition in the failure process occurs at this 
point. What is the nature of this transition? As we shall 
demonstrate, it is not a phase transition, but a crossover 
phenomenon. 

When localization sets in immediately in the break- 
down process, it is normally expected that the system is 
infinitely fragile in the limit of infinitely large system: As 
soon as a single fiber breaks at a given load, the entire 
system breaks down at that load [4| . As we will demon- 
strate, this is not the case here. 

We base our work on the discretized model for inter- 
facial failure proposed by Batrouni et al. [5]. A square 
array of L x L = TV linearly elastic fibers placed a dis- 
tance a apart connects a stiff half-space with a linearly 
elastic half space characterized by a Young modulus E 
and a Poisson ratio v (which we assume a typical value 
of 0.25 in the following). Each fiber, indexed by z, has an 
elastic constant k and fails irreversibly if it is elongated 
beyond an individual threshold value U. The threshold 
values are drawn from a uniform distribution on the unit 
interval. 

The separation of the two half spaces are controlled 
by displacing the hard medium by a distance D orthog- 
onal to the interface where the fibers sit. Fiber i then 
experiences a force 



fi = -k(ui - D) 



(1) 
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where ui is the local displacement of the softer half space 
at the position of fiber i. The forces from the fibers are 
transmitted through the softer elastic medium via the 
Green function [6|-[8j 
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12 \n(x,y) - rj(x',y')\ 



(3) 



fi — fj denotes the distance between fibers i and j at 
positions ?i and fj respectively. 

The Green function (|3|) is modified by the presence of 
boundaries due to the finite size (La) x (La) = Na 2 of 
the interface. We assume periodic boundary conditions 
and take into account the first reflected images. 

We note that if distances are measured in units of a, the 
Green function ([3]) is proportional to (Ea)~ x . Likewise, 
from Eq. (pQ), we see that the elastic constant of the fibers, 
k must be proportional to a 2 . Hence, if we change the 
linear size of the system, (La) — >> A(La), while keeping 
the discretization a fixed, we change only L —> XL in 
the model whereas we keep the parameters Ea and k 
fixed. On the other hand, if we change the discretization 
a — >• a/X while leaving the size of the system fixed, we 
change L — )• XL and the parameters (Ea) — )• (Ea)X and 
k -> k/X 2 . 

A given fiber breaks irreversibly (its elastic constant 
is set to zero) if stretched beyond a threshold value as- 
signed from a spatially uncorrelated probability distribu- 
tion. We choose the simplest: a uniform distribution. 
The model is quasi- static, and in lieu of time, we mea- 
sure the fraction of fibers that have broken, denoted by 
p. The load carried by the system is a(p) = ^2 i fi/N, 
and when a reaches its maximum, any extra load will 
result in a complete catastrophic failure. We denote this 
the critical load, cr Cl and the corresponding p c , the failure 
point. 

In the limit of (Ea/L) — >> oo, the model becomes iden- 
tical to the equal load sharing fiber bundle model (ELS) 
[US 03- On the other hand, for small values of (Ea/L), 
it does not approach any existing models. Models do ex- 
ist, e.g. the local load sharing fiber bundle model (LLS) 
PJJ], where the nearest surviving fibers absorbs the en- 
tire load that a fiber was carrying when failing. Another 
model, introduced by Hidalgo et al. [l2|, distributes the 
added load around a failed fiber as a power law in the 
distance from the failed fiber. In both models, there is 
no elastic response by the planes defining the interface. 

We have studied systems of size L = 256, L = 128, 
L = 64, L = 32 L = 16 and L = 8 with 10, 100, 1000, 
10000, 10000 and 10000 samples respectively. We explore 
a range of elastic constants e = (Ea/L) in the range 
e so ft < e < e st iff where e so ft = 2 -17 = 7.63 x 10 -6 and 
e stiff = 2 6 = 32. 

In order to visualize localization, we record the square 
distance between consecutively failing fibers, Ar 2 . If the 
positions of the failing are completely random, as is the 
case in the ELS fiber bundle model, the average distance 
is (Ar 2 ) 1 / 2 = L/y/6 « 0.408L. We show in Fig. □ a 
succession of histograms of (Ar 2 ) 1 / 2 . That is, we 

record Ar 2 (n) as a function of the number of failed 
fibers, n = pN — our "time" parameter. We then sum 
the number of times Ar(n) has had a particular value 
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FIG. 1. (Color online.) Histogram over the distance between 
consecutively failing fibers, (Ar 2 (n)) 1 / 2 , as function of the 
number of failed fibers n — pN. Darker colors correspond to 
higher density. The vertical bar indicates the failure point, 
p c N. The curve shows the running average (Ar 2 (n)) 1//2 . In 
all figures, L = 128. 100 simulations was used to construct the 
histogram. From top to bottom, e = 32, e = 2" 3 678 = 0.0781, 
e = 2~ 6 = 0.0156 and e = 2~ 17 = 7.63 • 10~ 6 . 



(Ar 2 ) 1 / 2 at n for several independent simulations, hence 
creating a histogram for each n. Darker colors signifies 
more hits at that value of (Ar) 1 / 2 . The curve shows the 
average value (Ar 2 (n)) 1 / 2 . With L = 128, we see that 
for the four different elastic constants e that we show, 
e = 32, e = 2 -3,678 (this value is chosen to make the 
figure comparable to the results in Batrouni et al. [Ej], 
where L = 128, E = 10, and a = 1, which gives e = 
0.0781 = 2 - 3 - 678 .) e = 2- 6 and e = 2~ 17 , (Ar 2 (n)) x / 2 
starts out being close to the ELF fiber bundle model value 
52.26. The vertical line in each figure shows the failure 
point n c = p c N. 

Rising the value if the elastic constant e above e st iff or 
below e so ft, will not lead to changes from the uppermost 
and lowermost panels in Fig. [TJ At the highest value 
of e, the system behaves as the ELS fiber bundle model 
throughout the entire breakdown process: The average 
distance between consecutively failing fibers, (Ar 2 (n)) x / 2 
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FIG. 2. (Color online.) Finite size analysis of failure point and 
critical loading. Both p c (main plot) and <j c (encapsulated) is 
plotted against l/log 10 [JV], for N = 8 2 = 256 to N = 256 2 = 
65536 and both have finite values for any N. The slope of the 
soft systems are a Pc =0.68 and a Gc — 0.16. 



remains constant throughout the process. The ELS fiber 
bundle model predicts p c = 1/2 so that n c = 8192. 

As the system gets softer, both (Ar 2 (n)) 1 / 2 and n c 
decrease. We see in three lower panels in Fig. [T] that 
there is an abrupt change in (Ar 2 (n)) x / 2 for some range 
of n values. This is localization. We also see that the 
failure point does not fall to zero as e is lowered. Even 
for the smallest value in Fig. [TJ e = e so ft = 2 -17 , n c is 
significantly different from zero. 

Fig. [2] shows the failure point Pc as a function of 
the inverse of the logarithm of total number of fibers, 
l/log 10 (iV). From this figure, we may extrapolate the 
value of p c in the limit of infinitely large system. We find 
p c (N — >• oo ) = p^ = 0.16 when e = e so ft and p^ = 0.5 
when e = e st iff. Likewise, we may extrapolate the critical 
load a c — see the insert in the figure. We find through 
extrapolation that <r£° = 0.17 for e = e so ft and <r£° =0.25 
for e = e s tiff, the value expected for the ELS fiber bundle 
model. 

It is a surprising result that neither p^ nor <r£° are zero 
for small value of e. The LLS fiber bundle model predicts 
that cr c - l/log 10 (7V) with a™ = 0, EH • Hidalgo et 
al. [HI present numerical evidence that their model also 
has <r£° = 0. Hence, in both of these models, = 0. 

Batrouni et al. |5j studied the structure of the clusters 
of failed fibers at p Cl claiming that at the failure point, 
they are distributed according to a power law with expo- 
nent -1.6. This would indicate a critical point at p = p c . 

Referring to Fig. [TJ we see that for e = e st iff = 32, the 
system behaves as the ELS fiber bundle model where the 
position of the fibers that fail bear no correlations among 
themselves. It is clear when observing the average dis- 
tance between consecutively failing fibers, (Ar(n) 2 ) 1 / 2 
which essentially remains close to the ELS fiber bundle, 
where (Ar 2 (n)) x / 2 = L/\/E. Hence, we expect that the 
clusters follow percolation theory [l5j. In Fig. [3] we show 
the density of of the largest cluster of failed fibers, s*, as 



FIG. 3. (Color online.) Density of the largest cluster of failed 
fibers, s*, as a function of the damage p for e = e s tiff and 
e = e so ft- Here L = 128. The vertical bar indicates the 
percolation point p = 0.59274. 
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FIG. 4. (Color online.) Fluctuations of p c , W c — 
yj{vl) ~ (Pc) 2 plotted against L for L = 8 to L = 256. The 
slope of the black line is —2/3. The stiff system follows the 
slope, while the soft systems deviates. 



a function of p for e = e st iff and e = e so ft- In the case of 
the soft system, we see that at p w 0.25, the largest clus- 
ter becomes visible and grows essentially linearly with p. 
This behavior is due to localization. When p approaches 
1, there are jumps in 5*, because of coalescence of clus- 
ters. On the other hand, when e = e st iff, we see behav- 
ior consistent with percolation theory. When p is in the 
vicinity of p = 0.59274, the site percolation threshold on 
the square lattice [16], s* shoots up and thereafter evolve 
linearly in p. 

The failure point at which the system fails catastroph- 
ically, p c , occurs long before the jump in 8* for e = e st iff. 
This is an indication that the system is not critical at the 
failure point. Schmittbuhl et al. P3 measured the fluc- 
tuations of the failure point as a function of the system 
size, finding Ap c ~ 1/L 0,65 . This is consistent with the 
GLS fiber bundle model. Daniels and Skyrme [18[ showed 
that the statistical distribution of the critical elongation 
in the GLS fiber bundle model has has the form 

p(u c )du c = N^ 3 f[CN^ 3 (u c - (u c ))]du c , (4) 

C is a constant only dependent on the threshold dis- 
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FIG. 5. (Color online.) The critical load, cr c , plotted against 
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for several system sizes. The insert shows the rescaled 



a* is plotted against log 10 (e/e*), demonstrating data collapse. 



tribution and u c is the critical elongation. This leads 
immediately to 

Ap c ~ (K - {u c )) 2 ) 1/2 ~ iV" 1 / 3 = L- 2 / 3 , (5) 

where we have used the assumption that the threshold 
distribution is uniform in the vicinity of (u c ) in relating p c 
to u c . From Fig. 2] we can see that the stiff system scales 
as L -2 / 3 , while the soft systems deviates. We conclude 
that we cannot detect spatial correlation in the failure 
process beyond an uncorrelated percolation process for 
e e st iff. 

We now proceed to study the onset of localization. Fig. 
[5] shows the critical load a c as a function of e for systems 
of size L = 32 to L = 256. As e decreases, we observe that 
a c drops from the ELS fiber bundle value, goes through 
a crossover and ends in a stable a c for each L in the soft 
regime. We define e* by setting 



°~c(e*) = - (cr c (e s tiff) + cr c (e S oft)) • 



(6) 



<W =cr c (e) 



(?) 



We show <r*(e) vs. log 10 (e/e*) in the insert in Fig. 
As the correction term —a ac /log 10 [N] — >• in the 
macroscopic limit, we know that the curve never grow 
past <r c (e*) = 0.25. The largest gradient in <j c e seem to 
converge around 



Aa c (e) 



Alog 10 [e/e*] 



0.03, 



(8) 



and the shape of the curve is kept. The onset of localiza- 
tion is not a phase transition, but a crossover: there are 
no divergences anywhere in the derivative of this curve. 

The following picture then emerges: For a given elas- 
tic constant, e, the breakdown process starts out as de- 
scribed by the ELS fiber bundle model. The spatial cor- 
relations between the failing fibers seems to be so weak 
that it can be described as an uncorrelated percolation 
process. If the elastic constant e is large enough, the sys- 
tem will undergo both the ELS fiber bundle model failure 
point and the percolation transition. Depending on the 
threshold distribution, the ordering of the two events, the 
ELS failure point and the percolation transition, may be 
reversed. With lower elastic constant, localization sets in 
and breaks off the ELS fiber bundle breakdown process. 
When localization sets in, all failure activity is then es- 
sentially limited to the rim of a growing cluster of failed 
fibers. The onset of localization is an instability, not a 
phase transition. 
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